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ESTIMATING  RUNUP  ON  BEACHES: 

A  REVIEW  OF  THE  STATE  OF  THE  ART 


PART  I:  INTRODUCTION 


Runup  is  the  movement  of  the  waterline  on  a  beach.  The  term 
"runup"  has  been  used  to  describe  the  wave  "uprush"  or  "swash"  on 
coastal  structures  for  at  least  several  decades.  It  was  first  the 
subject  of  an  entire  technical  paper  in  1953,  "Wave  Run-up  on 
Sloping  Structures,"  by  Granthem  (1953).  In  the  1950 's,  Thorndike 
Saville  and  others  at  the  US  Army  Coastal  Engineering  Research 
Center  conducted  an  extensive  series  of  tests  of  runup  on 
structures.  The  primary  focus  of  these  runup  studies  was  on  man¬ 
made  structures  such  as  dikes,  seawalls  and  revetments.  Coastal 
engineering  laboratories  around  the  world  have  continued  to  focus 
on  quantifying  wave  runup  on  structures. 

The  related  phenomenon  of  runup  on  beaches  has  not  received 
as  much  attention  in  the  engineering  and  scientific  literature. 
The  engineer  tasked  with  estimating  runup  on  a  beach  has  often 
adapted  some  of  the  available  methods  for  estimating  wave  runup  on 
structures.  This  usually  entails  assuming  the  beach  profile  is 
fixed  and  applying  Hunt's  (1959)  Equation  or  the  Shore  Protection 
Manual  (SPM)  (US  Army  1984)  curves  for  runup  on  structures. 
Saville 's  (1958)  method  for  composite  slope  structures  has  been 
used  in  attempts  to  account  for  the  non-uniform  slope  found  on 
sandy  beaches. 

Beaches  are  not  fixed  but  free  to  move  in  response  to  the 
incident  waves.  They  have  much  flatter  slopes  than  man-made 


3 


structures.  The  beach  slope  is  not  constant  across  the  profile 
and  offshore  sand  bars  which  effect  the  incoming  wave  energy  are 
often  present.  Thus,  the  attempts  to  apply  runup  methods  derived 
for  structures  to  beaches  have  been  educated  guesses  at  best. 


Importance  of  the  Problem 

Maintaining  existing  development  along  the  nation's  coastlines 
requires  improved  techniques  for  coastal  engineering.  The  extent 
of  runup  on  beaches  is  important  for  a  number  of  coastal 
engineering  and  coastal  management  issues.  The  limit  of  runup 
defines  the  zone  of  possible  wave  damage  and  thus  can  be  used  in 
the  determination  of  coastal  setback  lines,  the  design  of  the 
height  of  man-made  beaches  and  dunes,  models  of  beach  and  dune 
erosion  during  storms,  management  policy  related  to  the  creation 
of  dune  fields  on  developed  shorelines,  and  estimating  flooding  due 
to  overtopping  and  breaching  of  beaches  and  dunes . 

Beach  nourishment  is  the  most  widely  accepted  engineering 
technique  for  beach  management  along  much  of  the  nation's 
coastline.  Although  a  wide  sandy  beach  provides  protection  to 
landward  structures  during  major  storms,  quantifying  the  level  of 
protection  is  difficult.  During  storms,  sand  is  moved  alongshore 
and  offshore  with  the  net  result  often  being  erosion  of  the 
beachface  and,  if  present,  the  sand  dune.  Several  models  of  sand 
dune  erosion  in  response  to  storms  have  recently  been  developed. 
These  models  are  applied  in  the  design  of  beachfills  to  determine 
the  required  volume  of  a  "sacrificial"  sand  dune  to  survive  the 
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design  storm.  They  do  not  provide  information  on  how  high  the 
water  will  run  up  the  beach  or  how  high  the  dune  should  be  to 
prevent  overtopping. 

From  a  less  pragmatic  engineering  viewpoint,  runup  on  beaches 
has  been  investigated  to  describe  the  motion  of  runup  in  terms  of 
the  driving  force,  the  incident  waves.  Several  investigators  have 
described  runup  in  terms  of  its  frequency  spectra  and  discussed  the 
differences  between  the  runup  spectra  and  the  incident  wave 
spectra. 

The  most  important  dynamics  of  runup  on  beaches  can  be 
discussed  in  terms  of  the  effect  of  waves  on  a  popular  children's 
activity.  Sand  castles  built  below  the  high  tide  line  wash  away 
during  the  next  rising  tide.  The  individual  "wave”  (runup)  which 
does  the  most  damage  to  the  sand  castle  is  often  not  the  largest 
incident,  offshore  wave.  The  reasons  tor  this  lack  of  one-to-one 
correspondence  between  waves  and  runups  on  beaches  are  not  entirely 
clear.  Obviously,  the  interference  patterns  of  individual  runups 
on  the  beach  slope  are  part  of  the  explanation.  Backwash  down  the 
slope  from  the  previous  runup  can  decrease  and  even  eliminate  the 
runup  from  a  subsequent,  incoming  wave.  The  largest  individual 
runups  on  structures  often  occur  when  two  waves  combine  on  the 
slope  at  almost  the  same  time.  Thus  the  timing,  or  phase,  of  the 
individual  waves  on  the  beachface  is  critical  in  the  runup  process. 
This  interference  phenomenon  decreases  the  number  of  runups 
relative  to  the  number  of  incident  waves  and  changes  the  frequency 
spectra  of  the  runup  relative  to  that  of  the  waves. 

Surf  beat,  the  rise  and  fall  of  the  average  water  level  inside 
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the  surf  zone,  effects  extreme  runup  on  beaches. 


Surf  beat  has 


periods  between  30  seconds  and  several  minutes,  the  "infragravity” 
range  of  periods.  The  highest  runups  occur  when  the  inshore  water 
level  is  high  due  to  surf  beat  and  incoming  waves  can  propagate 
farther  inshore.  Surf  beat  has  been  neglected  in  the  estimation 
of  runup  on  structures,  however,  it  is  much  more  important  on 
beaches.  Many  nearshore  processes,  including  runup,  have  been 
explained  in  terms  of  the  energy  in  the  infragravity  range. 
However,  the  causes  and  effects  of  surf  beat  are  an  area  of  active 
research.  Although  the  sand  castle  problem  may  be  trivial,  the 
physical  processes  involved  are  the  important  processes  of  runup 
on  beaches . 

Purpose  of  this  Report 

The  specific  question  addressed  in  this  report  is,  given 
offshore  water  level  and  incident  wave  conditions  for  a  specific 
beach,  what  are  the  extreme  runup  levels? 

This  report  will: 

□  summarize  the  literature  on  runup  on  beaches, 

□  compare  the  different  methods  for  estimating  extreme  runup 
statistics  on  beaches, 

□  discuss  the  general  abilities  of  these  methods  to  provide 
engineering  guidance, 

□  provide  recommendations  for  engineering  applications  for 
beach  and  dune  design, 

□  discuss  areas  for  future  research  on  runup  on  beaches,  and. 


□  discuss  the  possibilities  of  future  applications  of  this 
research  to  engineering  questions  related  to  overtopping  of 
beaches  and  the  resulting  flooding  of  landward  areas. 

This  report  will  focus  on  the  engineering  question  of 
estimating  runup  extremes  on  a  given  beach  with  given  incident  wave 
conditions.  Thus,  this  report  will  take  a  rather  pragmatic 
engineering  approach  to  the  problem. 


Definitions 

The  vertical  distance  from  the  still-water  level  (SWL)  to  the 
intersection  of  the  water  surface  and  the  beach  slope  is  n  , 
shown  in  Figure  1.  The  time  history  of  n  ,  n(t),  can  be  plotted 
as  in  Figure  2.  The  mean,  n  ,  will  be  above  the  SWL  due  to  setup 
in  the  surf  zone.  One  definition  of  an  individual  runup  is  the 
vertical  distance,  Rp  ,  between  the  SWL  and  each  local  maximum,  or 
"peak,"  in  the  n(t)  record  (Figure  2a).  Another  definition  of 
an  individual  runup  is  the  maximum  vertical  distance,  Ru  ,  from 
the  SWL  between  successive  upcrossings  of  the  mean,  n  (Figure 
2b)  .  The  peak  definition  is  more  common  in  the  literature  on 
runup  on  beaches.  There  are  four  runups  by  the  peak  definition  but 
only  one  by  the  upcrossing  definition  for  the  n(t)  record  given 
in  Figure  2.  Individual  swash  height,  S  ,  is  defined  as  the 
vertical  minimum-to-maximum  distance  between  successive  upcrossings 
of  the  mean  in  the  r>(t)  record  (Figure  2c).  A]  1  notation  is 
summarized  in  Appendix  A. 
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PART  II:  LITERATURE  REVIEW 


This  literature  review  focuses  on  estimating  extreme  values 
of  runup  on  beaches.  Some  previous  influential  work  on  wave  runup 
on  structures  is  discussed  as  background  because  several  aspects 
are  applicable  to  runup  on  beaches.  Several  of  the  Dutch  papers 
present  models  of  equilibrium  profiles  of  gravel  beaches:  a 
component  of  these  models  is  the  runup.  The  literature  focusing 
on  the  spectra  of  runup  on  beaches  is  not  included  in  this  review. 
Oltman-Shay  and  Hathaway  (1989)  discuss  most  of  the  runup  spectra 
literature  in  their  "tutorial  on  infragravity  motions." 


Wave  Runup  on  Structures 

Hunt  (1959)  summarized  the  European  (Miche  1944,  1951; 
Iribarren  &  Nogales  1947)  and  American  (Granthem  1953,  Saville 
1957)  laboratory  data  available  on  wave  breaking,  reflection,  and 
runup  on  structures.  He  emphasized  an  equation  relating  runup 
height  to  wave  steepness  and  structure  slope  which  can  be  written 
as , 


R  /  H  =  (tan  a)  /  (H/LQ)0'5  (1) 

where  R  is  the  vertical  distance  of  runup  from  the  SWL,  H  is 
the  incident  monochromatic  wave  height,  LQ  is  the  deepwater 
wavelength,  and  a  is  the  angle  of  the  slope  of  the  structure. 

A  dimensional  form  of  Equation  (1)  is. 
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R  =  (HL0)05  tan  a  (la) 

Hunt  only  considered  monochromatic  waves,  i.e.  a  single  wave  train 
of  identical  height  and  period. 

The  general  form  of  Equation  (1), 

R  /  H  <x  (tan  a)  /  (H/L0)0'5  (2) 

is  referred  to  as  Hunt's  Equation. 

Battjes'  (1974)  doctoral  thesis  to  the  Delft  University  of 
Technology,  and  a  companion  paper  (Battjes  1972)  provide  both  an 
exhaustive  literature  review  of  runup  and  significant  original 
extensions  of  monochromatic-wave-based  methods  to  irregular  seas. 
He  emphasized  the  importance  of  the  right-hand  side  of  Hunt's 
Equation.  It  is  not  necessary  to  know  both  beach  slope  and  wave 
steepness  but  it  is  sufficient  to  know  only  the  combination, 

i  =  (tan  a)  /  (H/L0)0'5  (3) 

This  number  is  often  referred  to  in  the  literature  as  the  "surf 
parameter"  or  the  "Iribarren  number". 

Van  Oorschot  and  d'Angremond  (1968)  investigated  runup  on 
structures  due  to  incident  irregular  waves  at  the  Delft  Hydraulics 
Laboratory.  The  defined  runups,  Rz  ,  using  a  zero-upcrossing 
method  with  the  SWL  as  the  horizontal  baseline  for  determining 
upcrossing  instead  of  n  •  They  chose  the  runup  level  exceeded  by 
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2%  of  the  runups,  Rz2X  ,  as  a  measure  of  extreme  runup.  They 
related  Rz2t  to  a  form  of  Hunt's  Equation  modified  for  irregular 
waves, 

rz2*  7  Hs  =  C  (tan  a)  /  (H$/gTp2)0-5  (4) 

Hg  is  significant  wave  height  which  is  not  specifically  defined 
by  van  Oorschot  and  d ' Angremond  but  is  probably  the  average  of  the 
largest  one-third  waves  at  the  toe  of  the  slope,  Tp  is  the  period 
of  the  peak  of  the  frequency  spectra  of  the  waves,  g  is 
acceleration  due  to  gravity,  and  C  is  an  empirical  coefficient. 
They  found  C  to  be  a  function  of  the  width  of  the  wave  spectrum. 
Wider  incident  wave  spectra  generated  higher  extreme  runups.  They 
investigated  fixed,  impermeable  slopes  of  1:4  and  1:6. 

Kamphius  and  Mohammed  (1978)  investigated  runup  on  relatively 
steep,  smooth  structures  in  the  laboratory  at  Queen's  University 
in  Kingston,  Canada.  They  specifically  focused  on  wave  steepness 
and  beach  slope  combinations  which  would  yield  highly  reflective 
conditions.  Their  beach  slopes  ranged  from  1:1  to  1:3.  They 
found  Rza  to  be  about  2.4  times  the  average  runup.  This  is 
slightly  larger  than  would  be  predicted  by  an  assumption  of  a 
Rayleigh  distribution  of  runups.  They  noticed  a  slight  tendency 
for  the  ratio  to  increase  as  structure  slope  decreased. 

Ahrens  (1981)  combined  data  from  three  independent  studies  to 
describe  the  runup  statistics,  R^  ,  Rzs  ,  and  Rz  in  terms  of 
the  wave  steepness,  H$/gTp2  .  R^  ,  the  significant  runup,  is  the 
average  cf  the  largest  one-third  of  the  runups  and  R  is  the 


average  of  all  the  runups.  Ahrens  used  results  from  the  studies 
of  van  Oorschot  and  d'Angremond  (1968),  Kamphuis  and  Mohammed 
(1978),  and  Ahrens  (1979).  The  latter  study  investigated  runup 
due  to  irregular  waves  on  a  1:1.5  slope  in  the  laboratory  of  the 
US  Army  Coastal  Engineering  Research  Center  in  Fort  Belvoir, 
Virginia.  Ahrens  (1981)  presents  2nd-order  polynomial  regression 
curves  for  the  three  runup  statistics  for  slopes  ranging  from  1:1 
to  1:3.  He  also  presents  empirical  equations  in  the  form  of 
Hunt's  Equation  for  three  runup  statistics  for  a  1:4  slope, 

Rz2X  /Hs  =  1.61  E,  (5) 

Rzs/Hs  =  i-25  *  (M 

*2mm'**  -  °-84  *  (7) 

where  the  surf  parameter  is  defined  with  the  significant  wave 
height  at  the  toe  of  the  structure  and  the  deepwater  wavelength 
corresponding  to  the  period  of  peak  energy  density  of  the  wave 
spectrum. 

Mase  and  Iwagaki  (1984) 

Mase  and  Iwagaki  (1984)  conducted  laboratory  tests  of  runup 
due  to  irregular  waves  at  Kyoto  University.  They  used  impermeable 
slopes  of  1:5,  1:10,  1:20  and  1:30.  They  investigated  a 

number  of  different  aspects  of  runup  including  runup  spectra,  runup 
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distribution,  the  effect  of  wave  groupiness  on  runup,  and  the  ratio 
of  the  number  of  individual  runups  to  incident  waves. 

They  used  both  definitions  of  individual  runup  presented  in 
Figure  2.  The  ratio  of  the  number  of  runups  defined  by  the 
upcrossing  method,  Ru  ,  to  the  number  of  runups  defined  by  the 
peak  method,  Rp  ,  varied  from  0.4  to  0.9  depending  on  beach 
slope  (the  ratio  decreased  as  the  slope  decreased).  in  other 
words,  there  were  between  10%  and  150%  more  runups  using  the  local 
peak  definition.  R^  and  Rp  ^  were  about  10%  and  20%  less  than 
Rus  and  Ru  nean  '  respectively. 

Mase  and  Iwagaki  looked  at  the  decrease  in  the  number  of 
runups  compared  to  the  number  of  incident  waves.  The  ratio  of  the 
number  of  runups  to  the  number  of  waves  varies  from  0.2  to  0.9. 
Lower  ratios  occur  for  lower  £  ,  i.e.  for  flatter  beach  slopes. 

They  plot  the  logarithm  of  three  statist irg]  representations 
of  runup,  Rp  ^  ,  R^  ,  and  R|)ax  against  the  logarithm  of  wave 
steepness  for  each  slope.  The  data  fit  an  empirical  exponential 
relation  between  runup  and  wave  steepness  of  -0.37.  This  is 
different  than  the  exponent  of  -0.5  inherent  in  Hunt's  Equation. 
They  thus  present  two  different  empirical  formulations  for 
irregular  wave  runup  on  impermeable  slopes, 

R*/Hs  =  a  (tan  a)b  (HS/LQ)C  (8) 

R*/H$  =  d  [{tan  a)/  (Hs/Lo)05]e  (9) 
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where  R*  =  Rp  ^  R^,  Ru  „ean/  R^,  or  R(nax  .  The  empirical 
coefficients;  a  ,  b  ,  c  ,  d  ,  and  e  ;  for  these  equations  for 
both  definitions  of  individual  runup,  upcrossing  and  local  peak  are 
given  in  Table  1.  Although  neither  of  the  above  equations  clearly 
fit  the  data  better  than  the  other,  Mase  and  Iwagaki  preferred  the 
latter  equation.  Apparently  they  preferred  it  because  it  kept  the 
surf  parameter,  £  ,  intact  as  a  scaling  parameter. 

Mase  and  Iwagaki  found  that  the  highest  runup  often  occurs 
alone,  i.e.  not  followed  or  preceded  by  other  large  runups.  They 
attributed  this  to  the  interference  of  the  backwash  of  large  runups 
reducing  the  height  of  the  following  runup.  They  also  found  little 
dependence  of  the  runup  statistics  on  incident  wave  groupiness. 
They  attributed  this  to  the  fact  that  the  wave  groupiness  decreases 
due  to  the  breaking  process  before  runup. 


Mase  (1989) 

Mase  (1989)  adds  two  more  statistical  representations  of  runup 
to  his  empirical  model  from  Mase  and  Iwagaki  (1984).  He  includes 
coefficients  of  Equation  9  for  Rp2?  ,  the  level  exceeded  by  2%  of 
the  runups,  and  Rp(1/10)  ,  the  average  of  the  highest  1/10  of  the 
runups.  He  presents  these  coefficients  in  terms  of  the  local  peak 
method  of  defining  runup  and  does  not  give  the  corresponding 
coefficients  for  the  upcrossing  definition.  He  presents  plots  of 
predicted  versus  measured  Rp  ,  Rpa  ,  Rp(V10)  ,  RK  .  and  Rp  mx 
that  show  that  the  empirical  equations  fit  the  data  from  which  they 
were  calibrated  well.  The  correlation  between  predicted  and 


measured  values  is  0.98  to  0.99. 


Channell  et  al .  (1985) 

Channell  et  al.  (1985)  present  the  results  of  a  laboratory 
experiment  at  the  Hydraulics  Research  Station  in  Wallingford, 
England  on  runup  on  permeable,  shingle  (flat  gravel)  beaches  due 
to  monochromatic  waves.  The  mean  grain  size  for  the  "pea  shingle" 
used  as  the  beach  was  8.2  mm.  Initial  beach  slopes  ranged  from 
1:6  to  1:10  .  They  found  that  on  the  gravel  beaches,  runup  was 
about  one-third  of  that  predicted  by  Hunt's  Equation  for  smooth, 
impermeable  slopes.  They  found  this  result  both  for  waves  which 
were  too  small  to  significantly  redistribute  the  gravel  and  for 
larger  waves  which  did  move  the  gravel.  They  also  constrained  the 
gravel  in  wire  mesh  gabions  and  subjected  this  constrained  beach 
to  the  large  waves  and  found  runup  to  increase  significantly. 
Runup  on  the  constrained  beach  was  about  1.6  times  runup  on  the 
unconstrained  beach  for  the  same  wave  conditions. 


van  Hijum  and  Pilarczyk  (1982) 

Van  Hijum  and  Pilarczyk  (1982)  investigated  the  dynamic 
equilibrium  profile  of  gravel  beaches  in  the  Delft  Hydraulics 
Laboratory.  They  developed  a  model  for  the  shape  of  this  profile 
in  terms  of  wave  height  and  grain  size.  Their  model  also  includes 
the  angle  of  wave  approach,  wave  period,  and  initial  beach  slope. 
They  were  interested  in  runup  because  the  upper  limit  of  runup  is 
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the  limit  of  gravel  motion  and  the  limit  of  their  equilibrium 
profile.  For  monochromatic  waves  they  use  Hunt's  Equation.  They 
present  an  empirical  expression  for  the  coefficient  (=1  on  smooth 
impermeable  slopes,  and  less  on  real  beaches)  in  terms  of  the 
ratios  of  wavelength  to  grain  size  and  wave  height  to  grain  size. 
For  runup  due  to  irregular  waves  on  the  initial  beach  slope,  they 
present  another  empirical  equation.  Runup  is  related  to  grain 
size,  significant  wave  height  and  period,  and  initial  beach  slope. 
They  also  present  an  expression  for  wave  runup  on  the  equilibrium 
profile  which  is  a  function  of  wave  height  and  period  and  grain 
size  only  (oblique  angle  of  wave  approach  is  also  considered). 


van  der  Meer  (1988) 

Van  der  Meer  (1988),  in  his  doctoral  thesis  to  the  Delft 
University  of  Technology,  discusses  the  entire  spectrum  of  profile 
shapes  exposed  to  wave  attack.  At  one  end  of  the  spectrum  are 
static  structures  such  as  rubble  mound  breakwaters  with  large  armor 
units  that  do  not  move  due  to  wave  action.  At  the  other  end  are 
sandy  beaches  which  are  significantly  shaped  by  waves.  In  between 
are  dynamic  revetments  and  gravel  beaches.  The  larger  grain  sizes 
present  a  static  stability  problem  and  the  smaller  grain  sizes 
present  a  dynamic  stability  problem.  Van  der  Meer  presents  a  model 
of  the  dynamic  profile  of  the  beach  during  a  storm  which  is  a 
function  of  wave  height  and  period,  grain  size,  storm  duration, 
water  depth,  and  angle  of  oblique  wave  attack.  A  part  of  the  model 
is  an  expression  for  the  berm  crest  height  which  can  be  interpreted 
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as  an  estimate  of  extreme  runup  since  the  crest  forms  near  the 


limit  of  runup.  This  expression  is  given  as  a  function  of  wave 
height  and  period,  grain  size  and  number  of  waves.  The  number  of 
waves  appears  because  the  foreshore  slope  is  a  function  of  storm 
duration . 


Carlson  (1984) 

Carlson  (1984),  in  his  doctoral  thesis  to  the  University  of 
California,  investigated  runup  on  actual  beaches  on  San  Francisco 
Bay.  He  made  movies  of  the  incident  waves  and  the  resulting  runup 
at  two  different  locations.  He  used  standard  photogrammetric 
techniques  to  analyze  the  water  surface  elevation  and  the  location 
of  the  intersection  of  the  water  surface  and  the  beach  for  the 
moment  in  time  captured  by  each  individual  frame  of  photography. 
Ground  truth  was  provided  by  a  series  of  poles  extending  from  the 
beach  through  the  surf.  The  result  of  this  labor  intensive  method 
is  a  time-series  history  of  water  level  and  a  coincident  time- 
series  history  of  swash  location.  From  these  time-series,  Carlson 
computed  distributions,  frequency  spectra  of  runup  and  waves,  and 
coherence  between  the  two.  For  each  day,  a  time-series  of  waves 
and  coincident  swash  is  shown  graphically. 

Carlson  does  not  compute  either  wave  heights  or  runups  in  the 
standard  way.  He  never  considers  the  definition  (e.g.  zero- 
upcrossing)  of  an  individual  wave  or  runup  in  the  time-series  of 
water  elevation  and  swash  location.  He  analyzes  only  n(t)  ,  the 
vertical  distance  from  the  SWL  to  the  location  of  the  intersection 
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of  the  water  surface  and  the  beach.  Carlson  actually  measured  the 


horizontal  location  of  the  intersection,  which  could  be  called  nh , 
relative  to  an  arbitrary  baseline.  This  can  easily  be  converted 
to  the  corresponding  vertical  distance  from  the  SWL,  n  /  if  the 
beach  slope  is  known. 

The  wave  heights  measured  by  Carlson  were  very  small  for  field 
conditions.  Carlson  gives  the  root-mean-square  of  the  offshore 
water  surface  elevation  as  3 . 6  cm  and  10.0  cm  for  his  two  sets 
of  data.  Carlson  intentionally  did  his  experiments  on  days  of 
small  waves  to  avoid  any  low  frequency  energy  which  would  be 
present  in  rougher  seas.  The  surf  zone  was  very  narrow  (maximum 
width  of  roughly  3  m)  because  of  the  small  waves  and  the  steep 
beach  (1:9  and  1:12). 


Holman  (1986) 

Holman  (1986)  investigated  the  extreme  runup  statistics  on  an 
energetic  open-ocean  beach  at  Duck,  North  Carolina.  He  filmed  the 
swash  zone  during  three  weeks  in  October  1982  from  the  US  Army 
Engineer  Field  Research  Facility  pier.  The  individual  movie  frames 
were  digitized  and  analyzed.  Holman  designed  the  data  collection 
procedures  to  investigate  the  spectra  of  beach  runup.  He  was 
specifically  interested  in  the  energy  in  the  infragravity  range 
(this  focus  was  reported  in  Holman  &  Sallenger,  1985).  Holman  was 
later  approached  by  the  Coastal  Engineering  Research  Center  to 
reanalyze  his  data  focusing  on  parameterizing  extreme  runup. 
Holman  (1986)  presents  the  results  of  149  runup  time-series  with 
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incident  wave  heights  and  periods  ranging  from  0.4  to  4.0  m  and 
6  to  16  sec,  respectively.  Incident  wave  data  were  from  a 
waverider  buoy  in  20  m  depth  and  a  Baylor  gage  at  the  end  of  the 
pier  in  8m  depth. 

Holman  considered  four  runup  statistics;  Rpa;  na  ,  the  2% 
exceedance  value  for  the  time-series  of  n(t)  (sampling  rate  =  1 
Hz);  q  >  the  maximum  value  of  H(t)  ;  and  ,  the  2% 
exceedance  value  for  the  crest-to-trough  vertical  excursion  of 
individual  swashes  as  defined  by  the  upcrossing  method  (Figure  2c). 
Holman  found  that  each  of  these  four  parameters  could  be  explained 
in  terms  of  £  defined  as, 

K  -  a  /  (H^/L^05  (10) 

where  a  =  beach  slope  defined  as  the  average  slope  in  the  swash 
zone,  i.e.  the  beachface  or  foreshore  slope;  h^  =  spectral 
significant  wave  height  as  measured  by  the  wave  gage  in  about  8m 
depth;  and  LQ  =  linear  theory  deepwater  wavelength  corresponding 
to  the  period  of  peak  energy  density  of  the  wave  spectrum.  The 
tangent  is  dropped  from  beach  slope  term  in  Equation  10  since  for 
the  range  of  beach  slopes  found  on  natural  sandy  beaches,  tan  a 
-  a  (a  in  radians).  Holman  found  neither  wave  gage  was  superior 
in  explaining  runup  in  terms  of  Equation  10.  He  presents  plots 
and  regression  coefficients  of  each  of  his  four  runup  parameters 
versus  £  .  The  regression  coefficients  for  the  four  parameters 
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Incident 


with  and  without  setup  removed  are  given  in  Table  2.1 
wave  characteristics  are  from  the  Baylor  gage  in  8m  depth  for 
the  regression  coeffients  given  in  Table  2. 


Resio  (1987) 

Resio  (1987)  reanalyzes  Holman's  data  to  apply  an  extreme 
value  model  to  runup  on  beaches.  Resio  reexamines  Holman's  wave 
data  and  shows  that  the  two  sets  of  data,  in  8  m  and  in  20  m 
depth,  are  different  in  terms  of  wave  height.  During  some  data 
sets,  waves  increased  in  height  as  they  moved  into  shallower  water. 
During  other  data  sets  wave  heights  decreased  between  the  two 
gages.  Resio  claims  that  this  behavior  is  consistent  with  shallow 
water  wave  transformation  theory.  Essentially,  low  steepness  waves 
(swell)  waves  should  increase  due  to  shoaling  more  than  they 
decrease  due  to  bottom  friction  and  breaking.  For  high  steepness 
waves  (sea)  the  dissipative  processes  are  dominant.  Resio,  like 
Holman,  found  no  preference  for  either  of  the  two  gages  for 
explaining  runup  with  Hunt's  Equation. 

Resio  investigates  the  use  of  a  third  depth  for  wave  data 
collection,  1.5  m,  i.e.  well  inside  the  surf  zone.  He  finds  that 
these  waves  do  not  explain  the  Rp2J  statistic  in  terms  of  £  as 


’The  regression  coefficients  given  in  Holman's  (1986)  Table 
1  do  not  fit  the  data  plotted  on  the  corresponding  Holman  (1986) 
Figures  6d,  7a,  7b,  7c,  and  7d .  It  is  not  clear  whether  the 
coefficients  are  in  error  or  the  figures  are  plotted  incorrectly 
(Holman  1989,  personal  communication).  Table  2  of  the  present 
report  is  from  Holman's  (1986)  Table  la. 
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well  as  waves  at  the  two  deeper  water  gages.  Resio  suggests  that 
offshore  wave  data  explain  runup  better  than  inshore  wave  data. 
This  suggestion  is  not  supported  by  Resio 's  analysis  since  the 
shallow-water  gage  data  is  probably  not  representative  of  the  wave 
field  where  the  runup  was  measured.  The  gage  is  a  Baylor  gage 
attached  to  the  pier.  Large  scour  holes  along  the  pier  offshore 
of  the  gage  effect  the  waves  reaching  the  gege.  It  would  indeed 
be  paradoxical  if  Resio' s  suggestion  were  correct.  He  suggests 
that  it  is  preferable  to  go  farther  away  from  the  location  of  runup 
to  collect  incident  wave  conditions  to  explain  the  runup.  Unbiased 
data  such  as  could  be  provided  by  a  gage  immediately  offshore  of 
the  location  of  runup  measurements  are  necessary  to  make  such  a 
suggestion . 

Resio  applies  a  generalized  extreme  value  distribution  to 
runup  and  uses  Holman's  data  to  calibrate  the  model.  He  compares 
Carlson's  data  to  the  model.  Resio' s  approach  extends  Holman's 
approach  to  include  any  extreme  runup  value  of  interest.  Holman 
only  provides4  information  on  the  2%  statistic  and  the  extreme 
statistic  for  a  35-minute  duration.  Holman  did  not  attempt  to 
generalize  his  results  to  a  frequency  distribution  because  "it  was 
felt  that  the  small  number  of  data  points  was  insufficient  to 
select  the  appropriate  frequency  distribution"  (Holman  1986,  page 
532).  Resio  reanalyzed  Holman's  data  to  do  precisely  that.  Both 
used  the  same  basic  "model"  of  Hunt's  Equation  to  scale  the  runup 
phenomenon  with  the  exception  of  the  choice  of  wavelength  used  in 
the  surf  parameter. 

Resio' s  model  for  the  large  values  of  runup  (the  values  of 
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most  interest)  can  be  written  as: 


R'  =  1.25  -  1.05  (Tr-0.5)  _0-19  (11) 

or 

R'  =  1.25  -  1.05  [  (1/P0)-O.5]  -°-19  (12) 


where  Tr  =  the  return  period  of  the  runup  of  interest  in  number 
of  runups,  PQ  =  1/Tr  =  probability  of  exceedance,  R'  =  the 

8  r 

coefficient  of  a  form  of  Hunt's  Equation; 


(Rp*  -  Hl/H^  =  R '  E, 


(13) 


with 


(14) 


so 

(Rp‘  -  -05 

R'  -  - - -  -  (Rp  -  n)  a  (H^Lpi  °-5  (15) 

where  Rp*  =  is  the  runup  statistic  of  interest  as  defined  by  the 
local  peak  method,  =  incident  spectral  significant  wave  height, 
Lp  =  local  wavelength  corresponding  to  the  peak  period  of  the 
energy  density  spectrum.  The  location  of  the  incident  wave 
parameters  is  in  water  depth  of  8m.  Equations  11  and  12  are 


plotted  in  Figures  3a  and  3b,  respectively. 


Resio's  model  does  not  include  n ' s  effect  on  the  total  runup 
elevation.  Although  it  is  not  clearly  defined  in  Resio  (1987),  the 
dependent  variable  is  as  given  on  the  left-hand  side  of  Equation 
13  above. 


<RP'  -  d/h» 

The  n  has  been  removed  from  the  runup  elevation  (Resio  1989). 
This  is  a  serious  shortcoming  for  engineering  applications.  Holman 
was  unable  to  predict  r|  based  on  £  .  Although  n  is  a 
manifestation  of  wave  setup  in  the  surf,  it  is  not  the  usual 
engineering  definition  of  setup,  the  average  water  level  at  the 
intersection  of  the  SWL  and  the  beach. 

Two  examples  are  given  which  show  how  to  apply  Resio's  model 
with  an  assumed  n  •  The  examples  are  chosen  to  match  Holman’s 
conditions  for  comparison  of  the  two  models  in  the  next  chapter. 


**************************************************************** * 
EXAMPLE  PROBLEM  1: 

GIVEN:  incident  peak  wave  period,  Tp  =  10  sec 

incident  wave  height,  H^  =  2  m  in  8m  water  depth 

beachface  slope  =  0.10  (1:10) 

storm  wave  and  water  level  duration  =  35  min 

n  =  l  •  0  m 

FIND:  An  estimate  of  the  maximum  runup  elevation  during  the 

storm . 
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SOLUTION: 


1)  Calculate  the  probability  of  exceedance  of  interest. 

The  number  of  runups  during  the  storm  is 
N  =  ( 35  min)  (60  sec/min)  /  (10  sec/wave)  =  210  waves 

******************************************* 

*  Note:  Resio's  method  assumes  that  * 

*  (number  of  waves)  =  (number  of  runups)  * 

*  This  is  not  correct.  * 

******************************************* 

Therefore,  for  the  maximum  runup, 

PQ  =  1/N  =  1/210  =  0.048 

o 

2)  Calculate  the  Hunt's  Equation  coefficient  from  Eqn.  12  or  Figure 
3b, 

R'  =  0 . 872 

3)  Calculate  local  wavelength  (d  =  water  depth), 

L0  =  (g/2TT)  Tp2  =  (9.8  m/sec2)/2TT  (10  sec)2  =  155  m 
d/L0  =  (8  m)  /  (155  m)  =  0.0516 

from  lineai.  wave  tneory  tables.  Appendix  C  of  SPM, 
d/L  =  0.096 
so , 

Lp  =  d  /  (0.096)  =  (8  m)  /  (0.096)  -  83  m 

4)  Calculate  maximum  runup, 

R„x  ■  R’  °  <>V-)0-5  +  n 

=  (0.872)  (0.10)  [(2  m )  ( 8  3  m)]0'5  +  n  =  1-12  m  +  rj 

=  2.12  m 

5)  Calculate  elevation  of  maximum  runup, 

elevation  =  SWL  +  RmJIV  =  SWL  +  2.12 

max 

where  SWL  includes  the  tide  level  and  storm  surge. 
***************************************************************** 
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The  next  example  calculates  the  statistic  for  the  same 
given  conditions. 

***************************************************************** 
EXAMPLE  PROBLEM  #2: 

GIVEN:  The  same  storm,  wave  and  beach  conditions  in  Example 

Problem  #1. 

FIND:  The  elevation  exceeded  by  2%  of  the  runups  during  the 

storm  (runups  are  defined  by  local  peak  method). 

SOLUTION: 

1)  The  probability  of  exceedance  of  interest,  is  stated  as, 

P0  -  0.02  (2%) 

2 )  Calculate  R ' , 

R'  -  0.75 

3)  Calculate  Lp, 

Lp  =  83  m  (from  Ex.  Prob.  #1) 

4)  Calculate  runup, 

Rp2j  =  (0.75)  (0.10)  1(2  m)  ( 83  m)  ]CS  +  n  =  1-05  m  +  n 

=  2.05  m 

5)  Calculate  elevation  of  runup, 

elevation  =  SWL  +  2.05  m 

***************************************************************** 
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PART  III:  COMPARISON  OF  DIFFERENT  METHODOLOGIES 
FOR  ESTIMATING  RUNUP  ON  BEACHES 

Summary  of  Available  Methods 

The  methods  available  for  estimating  runup  on  beaches  are  very 
limited.  Holman's  method  and  Resio's  extension  of  Holman’s  work 
are  the  only  methods  based  on  actual  runup  measurements  on 
prototype  beaches.  Other  methods  are  based  on  laboratory  tests 
on  gravel  beaches,  impermeable  slopes,  or  revetments.  This  chapter 
will  focus  on  these  two  methods  and  the  laboratory  results  of  Mase. 
Mase  (1989)  and  Mase  and  Iwagaki  (1984  )  are  based  on  the  same 
investigation  and  will  be  referred  to  as  Mase's  work.  He  measured 
extreme  runups  due  to  irregular  waves  on  slopes  approaching  the 
flatness  of  natural  beaches  and  found  many  of  the  phenomena  of 
boech  runup. 

This  chapter  will  describe  how  to  use  Holman's  and  Resio's 
methods  for  engineering  applications.  Results  from  the  two  methods 
will  be  compared.  The  comparison  is  good  as  would  be  expected 
since  Resio  calibrated  his  method  with  a  subset  of  Holman's  data. 
Mase's  equations  compare  better  with  Holman's  data  than  previously 
reported.  The  problem  of  choosing  an  appropriate  beach  slope  for 
applying  Mase's  results  to  prototype  beaches  is  discussed. 
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Summary  of  Available  Data 

Independent  data  should  be  used  to  evaluate  the  available 
methods  of  estimating  beach  runup.  This  is  not  possible. 
Essentially,  no  beach  runup  data  in  either  raw  or  partially 
analyzed  form  were  available  for  this  investigation. 

Holman's  data  is  by  far  the  largest  set  of  beach  runup  data 
discussed  in  the  literature.  It  forms  the  basis  for  the  two  models 
of  beach  runup.  It  is  "not  available"  in  its  raw  form  of  digitized 
time-series  of  runups  and  collocated  beach  slope  and  incident  wave 
statistics  (Holman  1989,  personal  communication).  Holman's  summary 
statistics  of  his  149  runs  have  not  been  published.  Holman  (1986) 
only  published  dimensionless  plots  of  the  data  which  do  not  allow 
for  further  analysis  by  other  investigators.  A  summary  table  of 
Holman's  data  was  located  and  expanded  through  a  series  of 
telephone  conversations  with  Dr.  Robert  Holman  at  Oregon  State 
University,  Dr.  Donald  Resio  at  Offshore  and  Coastal  Technologies, 
Inc.  in  Vicksburg,  Mississippi,  Mr.  William  Birkemeier  at  the  Field 
Research  Facility,  and  Mr.  John  Ahrens  at  WES.  This  summary  table 
of  Holman's  149  data  points  is  given  in  Appendix  B.  It  includes 
the  location  and  time  of  each  data  run,  the  extreme  runup 
statistics,  computed  setup,  SWL,  beachface  slope,  and  incident  wave 
statistics  from  the  Baylor  gage  at  the  end  of  the  pier. 

Most  of  the  runup  data  were  collected  within  150  m  of  the 
pier  to  the  north  and  the  south.  The  data  may  have  been  biased 
because  the  wave  field  was  effected  by  the  perturbed  bathymetry 
around  the  pier  (Holman  and  Sallenger  1984).  Bathymetric  plots  of 
the  area  before  and  immediately  after  Holman's  experiment  are  shown 
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in  Figures  4  and  5.  The  pier  effects  on  bathymetry  extended  100 
m  to  400  m  to  either  side  of  the  pier.  Considering  the  natural 
directional  spread  of  energy  in  a  wave  field,  wave  energy  from  the 
direction  of  the  pier  must  have  been  modified  by  the  pier  effects 
for  all  of  Holman's  data. 

The  pier  appears  to  have  affected  the  runup  data  measured  by 
Holman.  Inspection  of  the  n  results  for  the  two  major  storms  in 
the  data  set  shows  this  effect.  These  two  storms  account  for  80 
of  Holman's  149  observations  and  for  all  of  his  larger  wave 
observations.  During  the  first  part  of  the  storms  (10-11  Oct  and 
24  Oct),  waves  and  winds  were  driving  the  longshore  current  to  the 
south  (Mason  et  al.  1984).  Later  in  both  storms  (12-13  Oct  and  25 
Oct),  waves  and  winds  came  around  and  longshore  current  reversed 
to  the  north.  n  was  higher  on  the  updrift  side  of  the  pier 
throughout  both  storms.  A  possible  explanation  is  the  effect  of 
the  bathymetry  near  the  pier  on  the  incident  wave  field. 

Carlson  published  plots  of  two  runup  time-series.  Resio  used 
these  to  check  the  parameter  estimation  of  his  generalized  extreme 
value  distribution  for  runup.  He  found  that  one  of  the  data  sets 
fit  his  distribution  and  the  other  data  set  may  have  been  biased 
since  the  incident  wave  conditions  were  measured  inside  the  surf 
zone.  Carlson's  data  set  does  not  reflect  open-ocean  beach 
conditions  since  he  measured  runup  due  to  very  small,  locally 
generated  waves  incident  on  an  unbarred  beach  profile.  Thus,  the 
agreement  with  Resio' s  model  is  surprising,  indeed  encouraging, 
considering  the  differences  in  the  overall  conditions. 

Guza  and  Thornton  (1982)  measured  runup  on  a  beach  using  gages 
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and  reported  on  the  runup  spectral  characteristics  but  not  the 
extreme  values.  It  is  not  known  whether  their  data  are  available 
in  a  raw  form  for  further  analysis.  More  runup  data  were  collected 
by  Holman  at  Duck,  North  Carolina  in  1985  and  1986  with  movie  and 
video  systems  (Mason  1987).  However,  analysis  of  the  data  is  not 
complete  (Holman  1989,  personal  communication ) . 


Mase's  Equation  Compared  to  Holman's  Data 

Both  Mase  (1989)  and  Holman  (1986)  state  that  Mase's  equation 
overpredicts  Holman's  data  by  100%.  Mase  says  that  his  equation 
is  roughly  an  upper  envelope  to  Holman's  data  and  that  dividing  his 
coefficient  by  two  provides  a  good  fit  through  the  data.  Mase 
attributes  the  difference  to  beach  permeability  but  does  not 
support  this  conclusion.  Neither  author  discusses  what  the 
appropriate  choice  of  beach  slope  for  applying  Mase's  results  to 
Duck,  North  Carolina.  Recall  that  Mase's  equation  is  based  on 
laboratory  data  on  constant  slopes  and  Holman's  data  is  prototype 
data . 

Figure  6  shows  a  beach  profile  at  Duck  before  and  during  the 
collection  of  Holman's  data.  The  depth  over  the  bar  is  1.1m  at 
mid-tide.  The  profile  is  reproduced  from  Holman  and  Sallenger 
(1985).  Holman's  wave  heights  were  0.4  m  <  H^  <  4.0  m.  Since 
common  values  for  the  ratio  of  spectral  significant  wave  height  to 
water  depth  at  breaking  are  around  0.5  ,  waves  were  breaking  on 
the  offshore  bar  for  most  of  Holman's  data.  The  surf  zone  during 
storms  included  all  of  the  profile  shown  in  Figure  6. 
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Superimposed  on  the  profile  in  Figure  6  are  slopes  of  a  =  0.3.00, 
0.050,  0.014  (1:10,  1:20,  and  1:70).  These  slopes  roughly 
correspond  to  the  beachface,  the  outer  slope  of  the  bar  and  the 
entire  surf  zone  during  major  storms. 

Figure  7  shows  how  Mase's  equation  for  R  compares  with 
Holman's  data  for  the  three  different  choices  of  beach  slope. 
Figure  7a  shows  the  comparison  of  Mase's  equation,  with  beach  slope 
of  m  =  0.1  (1:10),  with  the  envelope  of  Holman's  data  (from 
Holman's  1986  Figure  7c).  Figure  7b  shows  the  comparison  between 
Mase's  Equation  and  Holman's  data  using  a  beach  slope  defined  as 
m  =  0.05  (1:20).  The  agreement  is  good.  Figure  6c  shows  the 
comparison  between  Mase's  Equation  and  Holman's  data  for  a  beach 
slope  defined  as  m  =  0.0143  (1:70).  Mase's  Equation  underpredicts 
Holman's  data!  Thus,  depending  on  the  definition  of  beach  slope 
adopted  for  describing  the  Duck  profile,  Mase's  Equation  can 
overpredict,  correctly  predict,  or  underpredict  the  runup  measured 
by  Holman. 

Choosing  the  beach  slope  definition  to  characterize  the 
natural  beach  profile  of  Figure  6  is  a  problem.  Holman  used 
beachface  slope  to  define  his  surf  parameter.  His  data  collection 
and  analysis  technique  required  very  accurate  knowledge  of  beach 
face  slope  in  order  to  convert  from  horizontal  swash  data  to 
vertical  swash.  Thus,  he  had  good  beachface  slope  data. 

The  question  at  hand  when  trying  to  use  Mase's  results  is 
"what  single  beach  slope  should  be  used  to  correctly  model  the  Duck 
profile?"  By  inspection  of  Figure  6,  a  1:10  slope  should 
overestimate  runup  since  it  would  not  account  for  the  wave 
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transformations  including  breaking  and  dissipation  across  the  bar. 
A  1:70  slope  should  underestimate  runup  which  actually  occurs  on 
a  much  steeper  beachface.  Thus,  the  result  that  the  "best-fit" 
beach  slope  is  somewhere  between  these  two  extremes  should  not  be 
unexpected.  Note  that  the  slope  which  fits  the  data  best 
corresponds  to  the  slope  of  the  offshore  side  of  the  bar  in 
Figure  6 . 

The  beach  slope  definition  is  important  since  E,  has  been 
shown  to  be  a  very  good  scaling  parameter  for  a  number  of  surf 
phenomenon  (Battjes  1974).  However,  most  of  the  data  available 
during  the  development  of  the  surf  parameter  were  laboratory  data 
on  steep,  constant  slopes.  Since  £  is  a  scaling  parameter,  the 
experimentalist  is  free  to  define  the  exact  surf  parameter  in  terms 
of  the  wave  height,  length,  and  beach  slope  of  his  choice.  The 
surf  parameter  adapted  by  Holman  uses  estimates  of  the  three 
components  of  E,  from  three  different  locations  across  the  surf 
zone:  beach  slope  from  the  beachface,  wave  height  from  a  depth  of 
8  m  ,  and  a  representation  of  deepwater  wavelength.  Resio 
preferred  the  use  of  wave  height  and  wavelength  from  the  same 
location. 

Others  faced  with  engineering  applications  of  £-based 
methodologies  have  faced  the  same  problem  with  respect  to  beach 
slope.  Inman  and  Jenkins  (1989)  believe  that  swash  is  best  scaled 
with  beach  face  slope  and  surf  beat  is  best  scaled  with  the  overall 
beach  slope.  Dean  (1987),  in  the  keynote  address  to  the  Coastal 
Sediments  '87  Conference,  suggests  a  way  to  avoid  the  difficulties 
inherent  in  choosing  which  portion  of  the  beach  profile  is 
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appropriate  by  relating  surf  parameter  to  grain  size  and  incident 
waves . 


Resio's  Equation  Compared  to  Holman's  Equation 

Resio's  method  should  give  the  same  results  as  Holman's  method 
since  it  was  calibrated  with  a  subset  of  Holman's  data.  This 
section  will  show  that  the  two  methods  do  agree. 

The  results  of  Example  Problems  l  and  2  allow  a  rough 
comparison  of  the  two  methods  results.  The  35-minute  duration  is 
what  was  measured  by  Holman  and  the  beach  and  wave  parameters  fall 
within  the  ranges  observed  by  Holman.  Holman's  model  has  R' 
coefficients  of 

R'  =  0.65  (for  Rjgay)  /  and 
R'  =  0.53  (for  Rp2j )  . 

The  difference  between  Holman's  coefficients  and  those  from  Example 
Problems  1  and  2, 

R'  =0.87  (for  R|iax  ),  and 
R'  =  0.75  (for  Rpa)  , 

is  due  to  the  different  definitions  of  wavelength  employed.  Holman 
usea  deepwater  wavelength  and  Resio  used  local  wavelength. 
Conversion  can  be  approximated  since  wavelength  enters  Equation  14 
as  a  square  root.  For  this  problem  (d=8m). 
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L/Lq  =  83/155  =  0.54 

Converting  Holman's  coefficients  to  a  local  wavelength  definition 
of  psi  gives, 

R'  =  0.87  (0.54)0-5  =  0.64  (for  R^)  ,  and 

R’  =  0.75  (0.54)0  5  =  0.55  (for  R^} 

which  agree  well  with  those  from  Resio's  method.  Thus,  Resio's 
method  gives  consistent  results  with  Holman's  method. 

Resio's  method  allows  for  more  information  about  the  extremes 
to  be  estimated.  For  example.  Figure  3  can  be  used  to  calculate 
any  of  the  larger  runup  statistics.  Strictly  speaking,  Resio's 
model  should  be  limited  to  the  upper  10%  of  the  runups  since  that 
is  all  that  Resio  used  in  the  derivation.  However,  the  extreme 
value  distribution  used  by  Resio  is  related  to  the  Weibull 
distribution  and  others  have  found  that  a  Weibull  distribution  fits 
the  entire  runup  distribution.  Therefore,  the  lower  limit  of  the 
range  of  probabilities  for  which  Equation  12  is  applicable  may  be 
much  less . 


Summary  of  Comparisons 

The  few  appropriate  comparisons  between  available  runup 
estimation  methods  are  encouraging.  The  two  methods  based  on  field 
runup  data,  Holman's  and  Resio's,  agree  with  each  other  as  well 
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they  should  since  one  is  based  on  a  subset  of  the  other's  data. 
Resio  found  that  one  of  Carlson's  data  sets  fit  his  model.  The 
application  of  Mase's  laboratory  results  to  this  data  set  is 
complicated  by  the  problem  of  choosing  an  appropriate  beach  slope. 
For  a  beach  slope  intermediate  between  the  foreshore  slope  and  the 
overall  surf  zone  beach  slope,  Mase's  equation  predicts  Holman's 
data  well. 
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PART  IV:  ENGINEERING  APPLICATIONS  OF  BEACH  RUNUP  ESTIMATES 


The  development  of  models  to  describe  runup  on  beaches  may 
have  lagged  behind  the  development  of  models  to  describe  other 
coastal  phenomenon  because  of  the  difficulties  inherent  in  the 
phenomenon.  The  limit  of  runup  is  dependent  upon  a  number  of 
variables  including  the  water  level,  incident  wave  conditions, 
sediment  characteristics ,  and  the  beach  profile  which  is  constantly 
changing  in  response  to  the  waves  and  water  level.  This  problem 
of  the  changing  profile  or  cross-section  has  not  been  faced  by  the 
engineer  tasked  with  designing  a  traditional,  fixed  seawall. 
(However,  this  is  changing  with  the  recent  popularity  of  flexible 
revetments  designed  to  deform  without  loss  of  function  during 
storms.)  The  upper  limit  of  runup  would  seem  to  be  an  important 
component  of  coastal  engineering  and  coastal  zone  management 
decisions  since  this  limit  defines  the  zone  of  possible  wave 
damage.  This  chapter  will  discuss  several  of  the  areas  of  coastal 
engineering  where  estimates  of  runup  on  beaches  can  be  applied 
including: 

-  design  of  the  height  of  man-made  beaches  and  dunes, 

-  models  of  beach  and  dune  erosion  during  storms, 

-  management  policy  related  to  the  creation  of  dune  fields 
on  developed  shorelines,  and 

-  estimating  flooding  due  to  overtopping  and  breaching  of 
beaches  and  dunes . 

Engineers  may  realize  other  applications  if  given  a  reliable, 
practical  model . 
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Beach  and  Dune  Design 

Knowledge  of  the  limit  of  runup  could  be  applied  in  several 
ways  to  the  design  and  management  of  beaches  and  dunes.  Design 
engineers  tasked  with  prescribing  the  height  of  a  beach  berm  or 
dune  line  could  consider  both  the  design  water  level  and  the 
coincident  runup.  In  the  crudest  application,  the  height  of  the 
beach  or  dune  could  be  designed  as  for  a  fixed  structure,  i.e.  the 
required  height  to  prevent  overtopping  would  be  the  design  water 
level  with  storm  surge  plus  maximum  runup.  The  Federal  Emergency 
Management  Agency  (FEMA)  essentially  uses  this  concept  of  a  "fixed" 
dune  when  mapping  velocity  or  V-zones  for  flood  plain  mapping  in 
the  coastal  area.  The  delineation  of  these  v-zones  has  important 
local  economic  development  impacts. 

However,  since  storms  will  erode  the  beach  and  dune  and  a 
rational  design  approach  is  to  include  an  adequate  volume  of  sand 
in  the  dune  to  survive  the  design  storm  erosion  event.  Kriebel 
and  Dean  (1985)  and  Vellinga  (1953)  ha*,  e  developed  models  of  storm 
induced  dune  erosion  that  can  be  used  to  determine  the  volume  of 
the  required  "reservoir  of  sand."  The  State  of  Florida  uses  a  form 
of  Kriebel 's  model  in  determining  its  coastal  setback  line.  Such 
models  can  be  improved  by  including  runup  in  the  formulation.  In 
fact,  Kriebel  and  Dean  (1985)  conclude  that  the  first  improvement 
to  the  model  should  be  the  "inclusion  of  wave  runup"  on  the 
beachface . 

Knowledge  of  runup  on  beaches  could  be  applied  to  the  design 
and  management  of  grassed  dune  fields.  For  example,  the  minimum 
elevation  at  which  beach  grasses  and  small  dunes  can  be  maintained 
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on  the  beach  has  been  found  to  be  fairly  consistent  on  the  barrier 


island  of  Ocean  City,  New  Jersey  iWeggel,  1988).  This  minimum 
elevation  is  probably  related  to  some  recurrence  interval  of  runup. 
The  grasses  can  tolerate  only  a  certain  exposure  to  waves. 
Establishing  dunes  is  a  significant  coastal  zone  management  and 
policy  concern  on  many  developed  coastlines. 


"Overtopping"  of  Beaches  and  Dunes 

Overtopping  is  a  term  reserved  usually  for  the  water  passing 
over  a  fixed  seawall.  As  waves  break  on  an  overtopped  seawall, 
they  wash  up  (runup)  and  over  (overtopping)  the  structure. 
Knowledge  of  the  volume  of  water  overtopping  a  seawall  is  needed 
to  design  drainage  or  other  flood  prevention  systems  on  the 
landward  side.  Given  the  uncertainty  of  many  of  the  variables 
contributing  to  seawall  overtopping  our  ability  to  estimate 
volumetric  rates  of  overtopping  is  rather  poor.  Tn  a  review  of 
the  available  methodology,  Douglass  (1984)  concluded  that 
overtopping  estimates  were  probably  only  accurate  to  within  an 
order  of  magnitude.  Subsequent  laboratory  tests  with  i -regular 
waves  and  several  different  structure  cross-sections  at  WES  have 
greatly  improved  our  ability  to  estimate  overtopping  over  fixed 
structures.  With  structure-specific  laboratory  model  testing,  the 
accuracy  of  our  estimates  of  si te-speci f ic  overtopping  are  probably 
primarily  limited  by  the  inaccuracies  of  estimatinq  design  storm 
water  level,  the  single  most  important  overtopping  parameter. 
Prototype  overt oppi ng  data  are  needed  to  evaluate  the  importance 


of  scale  effects  in  the  laboratory  results. 

The  extension  of  the  concept  of  overtopping  of  structures  to 
overtopping  of  beaches  or  dunes  is  tenuous  and  should  be  approached 
with  many  caveats.  Essentially,  the  question  is  "what  happens  when 
the  runup  on  the  beach  exceeds  the  crest  of  the  sand  dune?"  Water 
from  the  individual  runups  flows  landward.  The  big  difference 
between  this  and  the  overtopping  of  structures  question  is  that 
sand  from  the  dune  will  be  carried  landward  with  the  water. 

Coastal  geologists  call  this  process  "overwash."  It  is  a 
mechanism  for  barrier  island  retreat  because  of  the  landward 
transport  of  sand.  Overwash  is  fairly  common  during  major  storms 
on  some  barrier  islands.  It  occurs  initially  at  low  spots  in  the 
dune  crest  elevation.  These  low  crest  areas  are  often  associated 
with  pedestrian  pathways  to  the  beach  where  the  vegetation  has  been 
killed  and  the  wind  has  created  a  "blowout"  in  the  dune  crest  line. 
The  width  of  these  "breaches"  in  the  dune  line  can  increase  during 
the  storm.  The  very  limited  quantitative  information  about 
overwash  is  primarily  limited  to  measurements  of  the  volume  of  sand 
transported.  Fisher  et  al.(1974)  and  Leatherman  (1977)  discuss 
some  simple  hydraulic  measurements  during  minor  overwash  events  on 
Assateaque  Island.  It  is  interesting  that  they  report  the  number 
of  runups  that  "surged"  through  the  breach  in  the  dune  averaged  out 
to  be  about  one  per  minute.  Although  they  don't  discuss  the  actual 
time-series,  the  timing  implies  that  these  "surges"  could  have  been 
driven  by  energy  in  the  infragravity  band. 

The  practical  problem  for  the  coastal  engineer  is  how  to 
estimate  flooding  on  the  landward  side  of  a  dune  when  the  dune  is 
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breached  or  overwashed.  There  exists  no  available  guidance  to 
address  this  problem.  Birkemeier  et  al .  (1987)  make  the  only 
feasible  recommendation  to  FEMA  concerning  overwash. 


"Because  the  flooding  will  be  caused  by  localized  dune 
breaching  and  overwash  (unless  the  surge  level  erodes  or 
overtops  the  dune  crest),  the  depth  of  flooding  is 
difficult  to  predict.  If  the  surge  is  of  short  duration, 
the  flooding  will  occur  as  weir  flow,  and  there  will  be 
a  drop  in  water  level  across  the  eroded  dune.  However, 
once  the  dune  overwashes  inland,  the  maximum  surge  level 
should  be  assumed  [as  the  interior  flood  level].  Since 
the  overwash  process  is  poorly  documented  and 
unquantified,  a  conservative  and  recommended  approach  is 
to  use  the  maximum  flood  level  if  the  dune  is  flooded." 


For  the  case  of  minor  overwash  events  it  is  tempting  to  try 
to  apply  overtopping  methodologies  developed  for  structures.  A 
model  for  this  condition  could  be  based  on  Battjes  (1974)  model 
for  overtopping  of  structures.  Battjes  assumed  that  each 
individual  runup  could  be  explained  in  terms  of  the  monochromatic 
form  of  Hunt's  Equation.  Battjes'  model  could  be  modified  by  using 
Resio's  (1987)  distribution  of  runup  on  beaches.  Input  to  the 
model  would  have  to  be  an  assumed  beach  profile  that  is  either 
fixed  at  a  low  elevation  or  changing  in  time.  It  could  represent 
a  typical  post-hurricane  profile  through  the  area  of  breaching. 
Design  SWL  must  be  less  than  the  highest  point  on  the  profile  for 
Hunt's  Equation  to  have  any  physical  meaning.  With  design  wave 
conditions,  Resio's  model  would  estimate  how  many  runups  would 
exceed  the  profile  crest.  The  results  of  Battjes  and  Roos  (1972) 
relating  the  runup  to  the  volume  of  water  on  a  slope  could  then  be 
used  to  estimate  the  overtopping  volume  rate  per  longshore  unit- 
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width  of  breach. 


Because  of  the  unverified  nature  of  such  an 


approach,  it  could  only  be  assumed  to  be  accurate  to  within  an 
order  of  magnitude. 


Research  Recommendations 

Research  on  runup  on  beaches  should  be  more  focused  on 
specific  engineering  applications .  This  can  be  accomplished 
without  sacrificing  the  general  applicability  of  the  research. 
Engineers  have  adapted  the  available  technology  to  their  specific 
engineering  problems.  It  is  telling  that  the  best  available  data 
set  on  beach  runup,  Holman  (1986),  was  not  collected  to  study  beach 
runup  for  engineering  purposes  but  in  order  to  study  the  swash 
spectra  with  a  mind  towards  identifying  the  importance  of  the 
infragravity  energy  in  the  surf  (Holman  and  Sallenger,  1985).  The 
engineering  applications  paper  derived  from  this  data  set,  Holman 
(1986),  was  a  reanalysis  of  the  data  that  was  commissioned  by  CERC. 
Although  the  application  of  available  technology  to  specific 
engineering  problems  is  necessary  and  productive,  engineering 
research  matures  to  studying  the  problem  of  interest.  The 
applications  for  beach  runup  models  discussed  above;  storm  erosion 
prediction,  beach  overwash  estimation,  dune  and  beach  design;  are 
important  and  current  questions  in  coastal  zone  management.  The 
use  of  beaches  and  dunes  for  storm  protection  is  not  new  and  will 
probably  play  a  more  major  role  in  the  nation's  response  to 
accelerated  sea  level  rise. 

Future  research  on  beach  runup  should  examine  whether  or  not 
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the  scaling  relationship  of  Hunt's  Equation  is  the  most  appropriate 
for  estimating  prototype  beach  runup.  Hunt  (1959)  developed  the 
equation  with  laboratory  monochromatic  wave  data.  His  data  showed 
a  clear  dependence  of  runup  on  both  the  wave  steepness  and 
structure  slope.  Wave  height,  wavelength  and  structure  slope  are 
clearly  defined  in  the  laboratory  with  monochromatic  waves. 
However,  the  definition  of  each  is  more  complex  on  prototype 
beaches . 

More  importantly,  the  importance  of  each  term  in  Hunt's 
Equation  on  beach  runup  data  has  not  been  clearly  examined.  Hunt 
(1959)  had  clear  evidence  that  runup  varied  roughly  linearly  with 
the  tangent  of  structure  slope  but  no  such  relationship  has  been 
shown  for  prototype  beaches.  Appendix  C  is  a  short  reanalysis  of 
Holman's  data  that  shows  that  runup  is  unrelated  to  beachface 
slope.  Only  one  site.  Duck,  North  Carolina,  is  included  in 
Holman's  data  set.  Data  from  other  sites  with  different  ranges  of 
slopes  must  be  considered  in  order  to  determine  whether  beachface 
slope  should  be  included  in  beach  runup  models. 

Mase's  results  support  the  idea  that  the  importance  of  beach 
slope  is  exaggerated  in  Hunt's  Equation.  Equation  8  separates  the 
slope  from  wave  steepness.  Mase  found  b=0.7  .  Thus,  Mase  found 
runup  varies  linearly  with  (tan  a)0,7  instead  of  linearly  with 
tan  a.  Mase's  b  is  fit  to  all  of  the  structure  slopes  he  tested 
(1:5,  1:10,  1:20,  1:30).  Examination  of  Mase  and  Iwagaki's  Figure 
10  shows  that  runup  varies  with  (tan  a)0J  for  the  two  flatter 
slopes.  In  summary,  slope's  influence  on  runup  decreases  for 
flatter  slopes. 
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Another  reason  why  the  scaling  relation  of  Hunt's  Equation 
should  be  questioned  is  that  the  dynamics  on  a  prototype  beach  are 
different  than  the  dynamics  on  simple  laboratory  slopes.  In 

particular,  overall  beach  slope  across  the  entire  surf  zone  width 
can  be  much  flatter  than  has  been  tested  in  laboratory  tests.  The 
presence  of  offshore  bars  is  very  common  on  prototype  beaches  and 
has  not  been  included  in  most  laboratory  tests.  Waves  break  on  the 
bar  and  propagate  across  the  bar  as  air-entraining,  turbulent  bores 
until  reaching  the  deeper  water  behind  the  trough  where  the  white- 
water  turbulence  ceases  and  one  or  more  waves  continue  propagating 
across  the  trough.  These  waves  which  have  reformed  in  the  trough 
are  the  waves  which  break  on  the  beachface  and  runup.  An 

individual  wave,  if  it  can  even  be  identified,  can  break  several 
times  over  a  distance  of  hundreds  of  meters  during  storm 
conditions.  The  final  breaking  wave  which  runups  ud  the  beach 
called  "shore  break,"  can  be  a  very  different  wave  than  the 
offshore  wave  both  as  an  individual  wave  and  in  terms  of  spectral 
description.  Similar  phenomena  have  been  described  on  beaches 
without  offshore  bars  (Guza  and  Thornton  1982). 

On  barred  profiles,  the  height  of  the  shore  break  is 
controlled  by  the  depth  on  the  bar.  The  influence  of  depth  across 
the  bar  on  runup  in  Holman's  data  is  investigated  in  Appendix  C. 
Dimensionless  depth  across  the  bar,  d/H,^  ,  explains  runup  as  well 
as  £  .  Since  depth  on  the  bar,  d  ,  was  not  measured  by  Holman, 
it  is  estimated  for  his  data  by  relating  SWL  to  the  profile  shown 
in  Figure  6.  That  such  a  crude  estimate  of  depth  across  the  bar 
explains  the  runup  as  well  as  E,  implies  that  further  work  on 
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finding  a  better  scaling  relationship  than  Holman's  application  of 
Hunt's  Equation  is  warranted. 

Research  on  beach  runup  should  consider  the  effects  and  causes 
of  energy  in  the  infragravity  portion  of  the  runup  spectra.  The 
rising  and  falling  of  the  water  level  in  the  surf  zone  allows 
incident  gravity  waves  to  propagate  farther  inshore  before  breaking 
and  beginning  the  runup  process.  Guza  and  Thornton  (1982)  found 
that  runup  spectra  is  saturated,  i.e.  that  the  amount  of  energy  in 
the  spectra  at  the  frequencies  of  the  incident  wind  waves  is 
independent  of  the  incident  wave  height,  while  the  energy  at  the 
infragravity  frequencies  continue  to  increase  with  increasing  wave 
height.  Mase  (1988)  found  the  same  results  in  his  laboratory  work 
and  showed  that  they  can  be  explained  entirely  in  terms  of  the 
interactions  of  the  waves  on  the  slope.  This  implies  that  a  runup 
model  separating  the  two  components,  call  them  swash  and  surf  beat, 
may  be  useful.  For  a  given  offshore  water  level  and  a  barred 
profile,  waves  are  depth-limited  across  the  bar  during  storms. 
Thus,  the  swash  component  of  runup  could  be  related  to  a  constant 
function  of  the  depth  across  the  bar.  This  would  essentially  be 
a  measure  of  the  wave  height  in  the  trough  which  could  then  be 
combined  with  the  beachface  slope  and  an  estimate  of  the  wave 
period  in  a  surf  parameter  which  should  be  related  to  the  swash 
height.  The  surf  beat  component  could  be  modeled  separately  and 
the  two  components  added  to  estimate  total  runup.  However,  the 
estimation  of  the  surf  beat  phenomena  on  beaches  is  still  an  area 
of  active  research. 

A  comprehensive  field  investigation  of  beach  runup  at  a 
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specific  site  could  be  modeled  after  different  portions  of  the 
recent  surf  zone  dynamics  experiments  at  Duck  (Mason  1987).  Waves 
and  water  levels  should  be  measured  offshore  and  at  several 
locations  across  the  surf  zone.  Specifically,  a  bottom-mounted 
wave  gage  in  about  8  m  depth  with  photopoles  or  bottom-mounted 
gages  on  the  bar  and  in  the  trough  are  a  minimum.  Photopoles  are 
painted  poles  driven  into  the  sand  to  extend  above  the  water 
surface  to  permit  movie  photography  or  video  recording  of  the  time- 
history  of  water  surface  elevation.  Runup  should  be  measured  with 
a  video  system  looking  in  the  alongshore  direction.  wind  and 
bottom  profile  data  should  also  be  available.  The  experiment 
should  be  conducted  away  from  the  local  bathymetric  effects  of 
piers,  groins,  and  other  coastal  structures.  All  of  the 
recommended  methodologies  and  the  required  analysis  techniques  are 
presently  used  by  CERC  so  no  new  equipment  or  methodologies  would 
have  to  be  developed.  Since  storms  are  of  particular  interest,  the 
problems  of  working  in  an  energetic  surf  zone  will  not  be  trivial. 

It  is  suggested  that  such  data  be  collected  at  more  than  one 
site.  In  particular.  Pacific  and  Gulf  coasts  should  be 
investigated  because  of  the  differences  in  storm  wave  climates. 

A  "package"  experiment  could  be  put  together  and  moved  to  several 
sites  around  the  country.  Experiments  should  be  scheduled  to  catch 
medium-sized  storm  events. 

Laboratory  tests  on  flat  slopes  and  non-uniform  slopes  are 
needed.  The  controlled  environment  of  the  laboratory  provides 
excellent  opportunities  to  investigate  specific  phenomenon. 
Laboratory  tests  are  one  of  the  most  powerful  coastal  engineering 
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tools  for  investigating  the  physics  of  complex  phenomena.  For 
example,  Mase  was  able  to  find  many  phenomena  in  his  laboratory 
flume  that  might  be  attributed  to  modes  of  propagation  of 
infragravity  energy  in  prototype  experiments. 
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PART  V:  SUMMARY 


phenomenon  of  runup  on  beaches  h^s  not  received  as  much 
attention  in  the  engineering  literature  as  the  phenomenon  of  runup 
on  structures.  In  the  absence  of  specific  guidance,  the  engineer 
tasked  with  estimating  runup  on  a  beach  has  often  adapted  some  of 
the  available  knowledge  of  runup  on  structures.  Holman's  method 
and  Resio's  extension  of  Holman's  work  are  the  only  methods  based 
on  actual  runup  measurements  on  prototype  beaches.  Holman  (1986) 
describes  a  large  data  set  on  beach  runup  and  provides  a  model  for 
estimating  extreme  values  of  runup.  Resio  (1987)  extended  Holman's 
model  by  describing  the  full  distribution  of  extreme  values. 

Holman's  data  is  published  for  the  first  time  in  Appendix  B 
of  this  report.  Limitations  of  the  data  are  that  they  are  only 
from  one  location  during  the  same  month  and  thev  mav  be  hiaRprt 
because  the  wave  field  incident  on  the  beach  may  have  been  affected 
by  the  perturbed  bathymetry  in  the  vicinity  of  the  pier. 

Both  Holman  and  Resio  used  the  same  basic  "model"  of  Hunt's 
Equation  to  scale  the  runup  phenomenon.  Resio  used  a  different 
choice  of  wavelength  in  the  surf  parameter.  Although  it  provides 
a  useful  method  for  the  generalization  of  Holman's  results,  Resio's 
model  is  of  limited  immediate  engineering  value  since  it  does  not 
yield  the  total  runup  elevation  and  it  does  not  account  for  the 
fact  that  the  number  of  runups  is  less  than  the  number  of  waves. 
Resio's  method  only  yields  the  variation  above  the  mean  runup 
elevation  (  n  )  which  is  difficult  to  estimate. 

Mase  (1988,  1989)  provides  the  most  significant  laboratory 
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contribution  to  understanding  the  dynamics  of  runup  on  baches. 
His  experiments  with  irregular  waves  on  sloping  structures  show 
many  of  the  same  properties  of  extreme  runup  statistics  and  runup 
spectra  found  on  natural  beaches. 

Future  research  on  beach  runup  should  closely  examine  the 
scaling  relationship  of  Hunt's  Equation  as  it  applies  to 
engineering  estimates  of  runup  on  beaches.  The  importance  of  and 
definition  of  each  term  in  Hunt's  Equation  has  not  been  addressed 
in  detail.  The  dynamics  on  a  prototype  beach  are  different  than 
those  in  laboratory  tests  with  monochromatic  waves  breaking  on 
planar  slopes.  In  particular,  overall  beach  slope  across  the 
entire  surf  zone  width  can  be  much  flatter  than  has  been  tested  in 
laboratory  tests.  The  effects  of  non-planar  beaches  on  wave 
breaking  has  not  been  adequately  considered  in  runup  models. 

Both  laboratory  and  field  experiments  are  needed.  Field 
experiments  should  be  modeled  after  the  recent  surf  zone  dynamics 
experiments  at  Duck,  North  Carolina.  However,  data  from  other 
locations  would  be  extremely  useful.  Focused  laboratory 
experiments  are  critical  in  that  they  can  provide  explanations  of 
the  field  results.  Mase's  work  is  an  excellent  example  of  how  the 
controlled  environment  of  the  laboratory  can  be  used  to  better 
understand  the  physics  of  beach  runup. 

Research  on  runup  on  beaches  will  be  used  in  the  solution  of 
a  number  of  coastal  engineering  and  coastal  management  problems. 
Since  the  limit  of  runup  defines  the  zone  of  possible  wave  damage, 
this  research  will  be  applied  to  the  determination  of  coastal 
setback  lines,  the  design  of  the  height  of  man-made  beaches  and 
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dunes,  models  of  beach  and  dune  erosion  during  storms,  management 
policy  related  to  the  creation  of  dune  fields  on  developed 
shorelines,  and  possibly  estimating  flooding  due  to  overtopping  and 
breaching  of  beaches  and  dunes . 
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Table  1.  Mase's  empirical  coefficients  for  Equations  (£) 
and  (9).  Adapted  from  Mase  and  Iwagaki  (1984) 

and  Mase  ( 1989 ) . 
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Table  2.  Holman's  (1986)  regression  coefficients  for  runup 

statistics  versus  £  as  given  by  the  incident  wave 
conditions  in  8  m  of  water. 
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Fi'  ’’re  1.  Definition  sketch  for  n  (t). 
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Figure  4.  Bathymetry  around  FRF  pier  at  Duck,  North  Carolina, 
24  August  1982  (from  Miller  et  al .  1986). 
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APPENDIX  A:  NOTATION 


a 
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c 

C 

d 
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g 

H 

H. 

Hs 
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N 
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P 
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R, 


max 


R„ 


R, 


p  mean 
Rps 


Rp(1/10) 


Empirical  coefficient  in  Equation  8 
Empirical  coefficient  in  Equation  8 
Empirical  coefficient  in  Equation  8 
Coefficient  in  Equation  4 
Water  depth,  m 

Empirical  coefficient  in  Equation  9 
Empirical  coefficient  in  Equation  9 
Acceleration  due  to  gravity,  m/sec2 
Wave  height,  m 

Spectral  significant  wave  height,  m 
Significant  wave  height,  m 
Local  wavelength,  m 
Deepwater  wavelength,  m 

Wavelength  from  linear  theory  corresponding  with  the 
peak  frequency  of  the  energy  density  spectrum,  m 

Number  of  waves  for  Resio's  model 

Subscript  signifying  local  peak  definition  of 
individual  runups 

Probability  of  exceedance  in  Resio's  model 

Vertical  runup  distance  due  to  monochromatic  waves,  m 

Maximum  runup  due  to  irregular  waves,  m 

Local  peak  definition  of  runup  (see  Figure  2),  m 

Mean  runup  with  local  peak  definition,  m 

Significant  runup,  average  of  the  highest  one-third 
runups  as  defined  by  the  local  peak  method,  m 

Average  of  the  highest  one-tenth  runups  as  defined  by 
the  local  peak  method,  m 


^7% 

Ru  man 
*us 

Ru2X 

Rz 

Rz  mean 

Rzs 

Rza 

R* 

R' 

S 

sa 

SWL 

t 

tan 

Tp 

Tr 

a 

n 
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Runup  level  exceeded  by  the  largest  2%  of  the  runups  as 
defined  by  the  local  peak  method,  m 

Upcrossing  definition  of  runup  (see  Figure  2),  m 

Mean  runup  with  upcrossing  definition,  m 

Significant  runup,  average  of  the  highest  one-third 
runups,  as  defined  by  the  upcrossing  method,  m 

Runup  level  exceeded  by  the  largest  2%  of  the  runups  as 
defined  by  the  upcrossing  method  m 

Upcrossing  definition  of  runup  relative  to  SWL,  m 

Mean  runup  with  SWL  upcrossing  definition,  m 

Significant  runup,  average  of  the  highest  one-third 
runups,  as  defined  by  the  SWL  upcrossing  method,  m 

Runup  level  exceeded  by  the  largest  2%  of  the  runups  as 
defined  by  the  SWL  upcrossing  method  m 

Generic  R  in  Equations  8  and  9 

Hunt's  Equation  coefficient  in  Resio's  model  (Equation 
15) 

Swash  height  (see  Figure  2),  m 

Swash  height  exceeded  by  the  largest  2%  of  swashes,  m 
Still  water  level 
Time,  sec 
Tangent 

Wave  period  corresponding  to  the  peak  frequency  of  the 
energy  density  spectrum,  sec 

Return  period  for  Resio's  model,  units  are  number  of 
waves 

Subscript  signifying  upcrossing  definition  of 
individual  runups 

Beach  or  structure  slope  in  radians  or  as  rise-over-run 

Vertical  distance  from  the  SWL  to  the  intersection  of 
the  water  surface  and  the  beach,  m 

Average  n  ,  m 


A2 


^max  Maximum  value  of  p  ,  m 

Level  exceeded  by  the  largest  2%  of  discrete  n(t) 
samples,  m 

H(t)  Time  history  of  n 

E  Surf  parameter 

Tt  Constant,  3.1416 
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APPENDIX  B:  HOLMAN'S  (1986)  DATA 


This  appendix  presents  the  summary  data  described  by  Holman 
(1986).  This  is  the  first  time  that  this  data  has  been 
published.  This  data  was  gathered  from  two  unpublished  sources. 
One,  a  summary  table  from  the  original  letter  report  from  Holman 
to  CERC.  Two,  a  set  of  two  large  binders  that  showed  some  of 
Holman's  analysis  results.  The  origin  and  a  brief  description  of 
each  column  is : 

1st  column:  Data  number.  Each  refers  to  a  movie  film  analyzed 
for  a  specific  transect  on  the  beach.  From  the  summary  table  and 
correlated  with  the  binders. 

2nd  column:  Date.  The  day  in  October  1983  when  the  movie  film  was 
shot  at  Duck,  NC.  From  the  binders. 

3rd  column:  Time.  The  time  of  the  17-minute  filming  (24-hour 
clock).  From  the  binders. 

4th  column:  Range.  The  distance  from  the  pier  of  the  runup 
analysis  in  meters.  "nH  and  "s"  refer  to  north  and  south  of  the 
pier.  Note  that  many  of  Holman's  149  data  points  are  from  the 
same  movie  film  at  different  locations  along  the  beach.  From  the 
binders  and  the  key  at  the  beginning  of  the  binders . 

5th  column:  setup.  This  is  actually  not  setup  but  r\  .  This 
column  of  data  is  from  the  summary  table.  Knowledge  of  the  tide 
gage  at  the  end  of  the  pier  must  have  been  used  to  calculate  this 
frcm  the  mean  of  the  n  record  which  was  referenced  to  an 
arbitrary  datum. 

6th  column:  Rs .  This  is  the  spectral  significant  runup  from  the 
summary  table.  From  the  summary  table. 

7th  column:  Hmo.  Spectral  significant  wave  height  from  the 
Baylor  gage  in  d=8m  at  the  end  of  the  pier.  From  the  summary 
table. 

8th  column:  T.  Tp  from  the  Baylor  gage.  From  the  summary  table. 

9th  column:  Tz .  The  mean  zero-uperossing  period  of  the  runups. 
From  the  summary  table. 

10th  column:  etamax.  rr  as  defined  in  Appendix  A.  From  the 
summary  table . 

11th  column:  eta2%.  as  defined  in  Appendix  A.  From  the 

summary  table  and  spot-checked  with  the  binders. 

12th  column:  Rp2%.  Ru2J  as  defined  in  Appendix  A.  From  the 
summary  table. 

13th  column:  S2%.  Sa  as  defined  in  Appendix  A.  From  the  summary 

B1 


table  and  spot-checked  from  the  binders . 

14th  column:  slope.  Beachface  slope.  From  the  summary  table. 

15th  column:  SWL.  Computed  with  the  n  from  the  summary  table  and 
the  first  moment  of  the  n  record  from  the  binders. 
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0.90 

9.5 

12.3 

1.26 

0.89 

1.03 

1.44 

0.14 

-0.30 

2 

4 

1555 

150s 

0.11 

1.47 

0.90 

9.5 

15.7 

1.15 

0.81 

1.02 

1.47 

0.12 

-0.30 

3 

5 

0836 

lOOn 

0.51 

2.10 

0.80 

11.0 

13.4 

2.08 

1.72 

1.94 

2.40 

0.14 

0.83 

4 

5 

1632 

250n 

1.29 

1.02 

0.80 

11.0 
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-0.35 
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5 
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200n 

-0.77 

2.94 

0.80 

11.0 

14.0 

1.21 

0.66 

1.04 

3.18 

0.16 

-0.35 

6 

5 

1632 

175n 

-0.42 

1.97 

0.80 

11.0 

12.2 

2.11 

0.69 

1.00 

2.41 

0.14 

-0.35 

7 

5 

1632 

150n 

0.15 

1.47 

0.80 

11.0 

13.9 

1.40 

1.02 

1.24 

1.72 

0.12 

-0.35 

8 

5 

1632 

I25n 

0.15 

1.65 

0.80 

11.0 

11.5 

1.30 

1.02 

1.19 

1.73 

0.14 

-0.35 

9 

5 

1632 

lOOn 

0.11 

1.56 

o.eo 

11.0 

13.0 

1.21 

0.93 

1.12 

1.58 

0.16 

-0.35 

10 

5 

1632 

lOOn 

0.05 

1.65 

0.80 

11.0 

12.1 

1.14 

0.8' 

1.02 

1.79 

0.16 

-0.35 

11 

5 

1632 

75n 

0.74 

1.08 

0.80 

11.0 

11.8 

1.57 

1.33 

1.41 

1.08 

0.15 

-0.35 

12 

6 

1024 

250n 

1.34 

1.48 

0.80 

11.0 

11.8 

2.38 

1.79 

1.81 

1.44 

0.11 

0.85 

13 

6 

1024 

200n 

1.02 

1.93 

0.80 

11.0 

11.5 

1.70 

1.87 

1.89 

1.95 

0.16 

0.85 

14 

6 

1024 

I75n 

0.92 

1.43 

0.80 

11.0 

12.4 

1.80 

1.64 

1.74 

1.48 

0.14 

0.85 

15 

6 

1024 

150n 

0.98 

1.41 

0.80 

11.0 

12.9 

1.86 

1.60 

1.71 

i.3i 

0.12 

C.85 

16 

6 

1024 

I25n 

0.97 

1.75 

0.80 

11.0 

14.4 

2.07 

1.72 

1.78 

1.70 

0.14 

0.85 

17 

6 

1024 

lOOn 

0.81 

2.33 

0.80 

11.0 

12.8 

2.19 

1.94 

2.06 

2.12 

0.16 

0.85 

18 

6 

1024 

75n 

0.84 

1.41 

0.80 

11.0 

11.7 

1.78 

1.41 

1.41 

1.42 

0.15 

0.85 

19 

6 
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50n 

0.62 

1.36 

0.80 

11.0 

11.8 

1.60 

1.39 

1.54 

1.40 

0.13 

0.85 

20 

8 
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0.54 

0.95 
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1.46 

1.15 
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8 

1308 
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1.17 
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1.34 
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22 
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23 

8 

1508 

150s 
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24 

9 
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0.49 

0.82 

0.42 
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0.98 
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0.86 
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0.77 

25 

9 
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0.07 
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0.42 

7.4 
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1.14 

0.65 

0.85 

1.04 

0.14 
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26 

10 
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1.01 

1.17 

2.17 

7.9 

12.9 

1.96 

1.67 
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0.11 

0.28 

27 

10 

0947 
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0.98 

1.48 

2.17 

7.9 

11.5 

2.60 

1.81 

2.00 
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0.16 

0.28 

26 

10 
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0.92 

0.94 
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7.1 
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1.87 

1.25 
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APPENDIX  C:  ANALYSIS  OF  HOLMAN'S  DATA 

This  Appendix  is  a  quick- and-d i rty  dimensional  analysis  of 
Holman's  data  which: 

1.  shows  there  is  little  influence  of  slope  on  runup. 

2.  shows  a  crude  estimate  of  dimensionless  depth  across  the  bar 
explains  the  data  as  well  as  £  does. 

3.  presents  an  empirically  modified  Hunt's  Equation  model. 


All  of  the  following  is  based  on  Holman's  statistic  which 
is  referenced  to  the  SWL.  Wave  information  appears  to  be  from  the 
Baylor  gage  at  the  end  of  the  pier  in  8m  depth. 

Figure  Cl  shows  the  dependence  of  runup  on  beachface  slope  is 
very  weak.  This  implies  that  the  wave  steepness  term  is 
responsible  for  the  ability  of  Hunt's  Equation  to  explain  runup. 
Figure  C2  shows  the  influence  of  Hunt's  wave  steepness  term  on 
runup. 

Depth  across  the  bar  was  estimated  by  considering  the  bar  to 
be  l.lm  below  the  SWL  for  all  cases  and  then  adding  the  n  for  that 
observation.  The  l.lm  is  probably  an  underestimate  of  the  depth 
of  the  bar  but  the  addition  of  "setup"  measured  on  the  beachface 
is  an  overestimate.  Figure  C3  shows  the  dependence  of  runup  on 
depth  across  the  bar.  The  relationship  looks  as  strong  as  that 
found  by  Holman  for  £  (Holman's  Figure  6b). 

Figure  Ci  combines  the  twu  dimensionless  terms,  me  equation 
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:  s  fit  to  this  data.  Thus,  a  resulting  model  of  runup  is 
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This  is  plotted  on  Figure  C5  and  the  following  equation  ran  be  used 
to  fit  Holman's  data. 


ir1  =  °-4  +  °*lum 
mo 


(Eg.  4) 


Cl 


where  a  modified  surf  parameter  is, 
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The  ability  of  the  above  equation  to  describe  Holman's  data 
is  shown  in  Figure  C6 .  Figure  C7  shows  how  well  Holman's  model 
fits  his  own  data.  The  present  mcdel  fits  better  but  has  three 
fudge  factors  (empirical)  instead  of  two,  so  it  should. 
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